Control of atomic decay rates via manipulation of reservoir mode frequencies 
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We analyse the problem of a two-level atom interacting with a time-dependent dissipative envi- 
ronment modelled by a bath of reservoir modes. In the model of this paper the principal features 
of the reservoir structure remain constant in time, but the microscopic structure does not. In the 
context of an atom in a leaky cavity this corresponds to a fixed cavity and a time-dependent external 
bath. In this situation we show that by chirping the reservoir modes sufficiently fast it is possible 
to inhibit, or dramatically enhance the decay of the atomic system, even though the gross reservoir 
structure is fixed. Thus it is possible to extract energy from a cavity-atom system faster than the 
empty cavity rate. Similar, but less dramatic effects are possible for moderate chirps where partial 
trapping of atomic population is also possible. 

PACS numbers: 42.50.Ct, 03.65.Yz 



I. INTRODUCTION 

The decay of an excited atomic state is manifestly af- 
fected by its environment. This has been evident since 
studies began on the effect of cavities and conducting 
plates on spontaneous emission [l|, [3, [3, 0] ■ The bound- 
ary conditions of nearby surfaces affect the electromag- 
netic mode density and, in the simplest treatment, the 
application of Fermi's golden rule shows that the spon- 
taneous emission rate is modified. However, with the 
development of high-Q cavities in the microwave [s*] and 
optical regime [6] it was possible to have an oscillatory 
exchange of energy between a single atom and a single 
cavity mode of electromagnetic (EM) radiation. In this 
regime a single quantum of EM energy may be involved 
and the atom may become entangled with the cavity. 
It is often desirable to control quantum systems of this 
kind, for example, to try and make a specific state of the 
cavity field (see e.g. Q), or to make entanglement, or in 
more complex systems, to make photons on demand Q. 
However, the effect of the environment is always there as 
decoherence, or dissipation, and it almost always acts to 
degrade the quantum phenomena in which we are inter- 
ested. 

In this paper we will look at a model example of what 
we might do in the environment to try and regain control. 
We will take the simplest quantum optical system of a 
two-level atom coupled to a leaky cavity. In our idealised 
system the cavity leaks through lossless dielectric mir- 
rors, although there well may be other systems that fit 
the model such as coupled microwave cavities ^ or engi- 
neered reservoirs for ion traps . For the optical cavity 
system, there are a number of well established types of 
quantum optical environment. For example, the thermal 
bath contains energy which ultimately leads to thermali- 
sation of the atom and cavity. Rigged reservoirs ^lT| have 
bath states with coherence which add gain and noise to 
the field state. The broad-band squeezed vacuum [l^l 
has correlations between different frequency bath states 
which in general has both the properties of a thermal 
bath and the presence of coherence. All of these envi- 
ronments would eventually excite an initially unexcited 



atom in the cavity. The environment we will look at is 
one in which the bath, or reservoir modes, change their 
frequency as a function of time. This kind of environment 
would not excite an initially unexcited atom, although 
the dynamics of an initially excited atom would be af- 
fected. In general, we would expect that reservoir modes 
do have time-dependent properties, since they are them- 
selves embedded in an uncontrolled environment. How- 
ever, here we will suppose that we can at least control 
their frequency for a short period of time. The reservoir 
will be fully specified if the frequency and the strength of 
the coupling to the atom are known for each mode: this 
amounts to specifying the reservoir mode structure. 

If we modulate all the reservoir mode frequencies by 
the same amount, keeping the coupling of each mode to 
the atom the same, the overall effect is the same as mov- 
ing one of the cavity mirrors, which has been studied for 
the case of an oscillating mirror • However, this case is 
also equivalent to a time-dependent detuning of the reser- 
voir from the atomic transition frequency. This situation 
has been discussed in Ref. where it was found that it 
was possible to inhibit the decay of the atom for partic- 
ular modulation frequencies. Modulation of the phase of 
the atom-cavity coupling constant was also analysed by 
Agarwal [iBj who showed that, for a low-Q cavity situa- 
tion, the decoherence could be reduced with fast modu- 
lation. Temporal variation of both the phase and weight 
of the coupling strengths between an atom and a contin- 
uum, with variation of the atomic detuning was studied 
in Ref. In this case the bath modes had to be mod- 
ulated identically and the time-dependence had to be a 
weak perturbation. However, it was shown that both en- 
hancement and suppression of the decay rate relative to 
the golden rule limit is possible. 

In this paper we aim to investigate a new and sub- 
tle case, where the mode frequencies and coupling con- 
stants are manipulated, but in such a way as to leave all 
macroscopic properties of the reservoir unaffected. The 
unaffected properties will be the overall spectral width 
and resonance frequency of the reservoir structure. To 
achieve this the bath modes both increase their frequency 
and change their coupling strength as shown in figure [T] 



2 



9 




g 




FIG. 1: Reservoir mode couplings gk{t) (see equation [1} at 
two diflerent times t\ and t2. The envelope of the couplings 
remains the same as the reservoir mode frequencies vary in 
time. Thus the coupling of the atom to an individual mode 
changes in time. The thick line indicates the same mode k at 
two different times, t\ and t2- Since at time i2, this mode has 
moved to a different frequency under the static envelope, the 
coupling has changed accordingly. 



Passing through the cavity resonance from below, the 
coupling strength of a mode to the atom first rises and 
then falls as the frequency is increased. The physical re- 
alisation could be a double cavity: one cavity couples to 
an atom and leaks to a second cavity which can have its 
end wall moved. We will discuss this realisation further 
in section Ivl 

Without further analysis this situation would present 
us with a dilemma. The macroscopic structure of the 
atom-reservoir interaction is unchanged in time, but the 
individual bath modes only interact with the atom for a 
finite time. It may not be immediately clear that this 
situation can dramatically affect the atomic decay: at 
any point in time the atom has the same coupling to 
some bath modes. However, because the bath modes 
change, the history of the atom-bath interaction affects 
the time development of their population. We will find 
that in this case both enhanced and inhibited decay are 
possible. 

The remainder of this paper is organised as follows. In 
section In] we describe our model for the problem and de- 
rive the equations governing the dynamics of the atomic 
state populations. Section Hill reviews the standard prob- 
lem of a two-level atom interacting with a high-Q cavity 
in the strong and weak coupling regimes. However, the 
section pays particular attention to the evolution of the 



quantum state of the reservoir degrees of freedom, which 
are often traced over in similar calculations. Section IIVI 
summarises the main results of the investigation in the 
regimes of strong and weak coupling. Finally, in section 
Ivl we estimate the size of the effects that might be ob- 
served and summarise the findings of the paper. 

II. DYNAMIC RESERVOIR MODEL 

We consider a two-level atom, with states denoted |0} 
and coupled to a reservoir of electromagnetic radi- 
ation modes at zero temperature. The reservoir modes, 
labelled by an index fc, have time- dependent frequencies 
ujk{t), and raising and lowering operators bl and bk- The 
Hamiltonian for the evolution of the composite atom and 
reservoir system is given in the rotating wave approxima- 
tion by {h — 1) 

H{t) = ujo^+a- + J2 ^k{t) {blbk + I/2) 

k 

+ J29kit){^'bl + bka+), (1) 

k 

where the atomic transition frequency is ojq, and = 
1)(0| and a~ — |0)(1| are the atomic raising and lower- 
ing operators respectively. The electric dipole coupling of 
the fc*'' mode of the electromagnetic radiation field to the 
atomic transition |0) ^ |1) is denoted by gk{t)- With- 
out loss of generality, the gk{t) are chosen to be real. 
For simplicity, we will later choose the spectral profile of 
the coupling constants to be peaked at the frequency of 
the atomic transition. While the reservoir modes are re- 
distributed, we will suppose that this envelope remains 
fixed. Therefore, the magnitude of the coupling constants 
gk{t) must vary in time to ensure this (see figure [T]). We 
will also assume that whilst the frequency change of the 
reservoir modes may be very rapid from the point of view 
of the atom, it will not be so fast that it could create pho- 
tons in the cavity by itself. In this sense the cavity field 
will adiabatically follow the motion of its end-mirror (see, 
e.g., Refs. ^m)- 

To develop equations for the state vector we work in 
an interaction picture and so we may write the time- 
evolution operator as 

U{t,0) = exp(^-i^ ^c^fe(T) (6^6fc + l/2)dr^ 

xexp(^-ioJota+a-^Ui{t,0). (2) 

The effect of the atom-reservoir interaction is contained 
in the term Ui{t, 0), which satisfies the Schrodinger equa- 
tion 

^^^4l^=mm,0), (3) 
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and the interaction picture Hamiltonian Hi(t) is given 
by 



gk{t)a cxp ( i j [^^(t) - cjojrfr 



+ gk(t)&^bk exp ( J [ajkir') - uJo]dT' 



(4) 



The state vector is expanded in terms of the eigenstates 
of the excitation number operator, N = &^a~ ^I^fc- 
Since N commutes with the Hamiltonian the num- 
ber of quanta is a constant of the motion. For initial 
states we assume the system is prepared with the atom 
in its excited state, and the reservoir is initially empty, 
i.e. IV'(O)} = |1) O I . . . . . .). Then the state vector at 
time t may be written as 



=Ca(t)|l)® |...0...) 



(5) 



where Ca(t) is the amplitude for the atom to be found in 
its excited state, and Ck{t) is the amplitude for the fc*'' 
mode of the reservoir to contain one excitation. Elimina- 
tion of Hi{t) and Ui{t) from equations ([3]), @ and ^ 
leads to the following coupled differential equations for 
the amplitudes Ca{t) and Ck{t): 



.dCa{t) 



.dckjt) 
' dt 



-'^ gk (t) exp (^-i J [ivkir) - ujo]dT 



Ck{t) 



(6) 



=gk (t) exp J [cdkir) - ujo]dT^ Ca(t). (7) 



We can now eliminate the amplitudes Ck (t) by integrating 
equation ([7]) and substituting the result into equation ([6]). 
Thus we obtain the integro-differential equation for Co(t), 



dt 



dt'Y, 9k{t)gk{t') 



X exp y-i j [uJk{T) - i^vi\dTj 

X exp i / [wfc(r') - wo]dT' Ca{t'). 



(8) 



For much of the following analysis we will work in the 
limit of a continuum of reservoir modes so that the sum 
in ([8]) is replaced by an integral over a continuous variable 
w. Thus we make the identifications 



LOk{t) 

9k{t) 



^g{uj,t) 



(9) 
(10) 



where uj refers to the initial frequencies of the field 
modes. The function ujji{uj,t) represents, at time t, the 
frequency of a bath mode which had frequency ujatt — O. 
Hence, we will have u)r{uj, 0) = uj. We also introduce the 
density of states p{uj) at time t = and then equation 
dH]) is replaced by 



dcgjt) 
dt 



dt' I duj p{uj)g{uJ,t)g{Lj,t') 



X exp ^-i y [ujii{uj,T) - ujo]dT^ (H) 
X exp i / [ujr{uj, t') - ujq] dr' Ca{t')- 



In addition, in equation (jlip . we have extended the lower 
frequency limit from to —oo. However, for optical 
transitions, (wq ^ lO^^s"^), the effects of this approx- 
imation are similar to others already made (for example 
the rotating-wave approximation, or the assumption of a 
two- level atom), and may be neglected. 

Later, we will also examine an idealised time- 
dependent bath spectrum, which is in fact just the fre- 
quency dependent population of the reservoir modes. 
We define this spectrum by considering those modes /ca 
(/cA G {fc}) which have frequencies that lie within Aw of 
a given frequency w at a time t. That is, we let 



S{io,t) 



(12) 



fcA 



In the continuum limit the number of modes /ca in the 
sum is approximately Aw/5(a;), and since the Ck are ex- 
pected to vary smoothly with k in this limit, we finally 
let 



(13) 



S{uj,t) — > p{uj)\ck^{t)\ 



which applies to a representative /ca- Equation ()13|1 will 
serve as an operational definition of the bath spectrum. 

We will now specify a model reservoir structure to be 
used in equation (jlip . For a simple cavity model, we can 
assume that the single time reservoir structure function 
p(ijj)\g{LO,t)\'^ is a Lorentzian, with width 7, centred on 
the atomic transition frequency wqi i-e. 



p{uj)\g{uj,t)\ 



7^ + {i^B{i^,t) - Wo)^ 



(14) 



This describes the time-dependent coupling of a bath 
mode, initially at to. However, it also describes the in- 
stantaneous coupling to the bath of modes at time t. 
Since the structure has no time dependence, apart from 
the time dependence of w^, the basic reservoir structure 
is constant in time, in this model. This is also clear from 
the weight D of the Lorentzian, which is defined through 
the relation 



p{u;)\g{Lo,t)\''dLO^D''. 



(15) 
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We can see that this integral is independent of time 
for the chosen form p^ . In the static case, where 
p{uj)\g(Lu,t)\'^ is not a function of time, the Lorentzian 
(1141) is a common choice for the reservoir structure 
[l9l [20| and the resultin g dy namics have been well ex- 
plored (see, for example, HI, In the weak coupHng 
limit, the Lorentzian reservoir structure ensures exponen- 
tial decay of atomic population. Of course, equation (fTTj) 
does not feature p{u])\g{uj,t)\'^ , but the two-time prod- 
uct p(uj)g(uj,t)g{uj,t'). However, as already mentioned, 
g{uj, t) may be chosen to be real. Therefore the two-time 
product follows immediately as 



we now substitute (I17p and ([T6|l into pTjl we find 



dcait) _ D^-f /•* ^^cxp [-Z (lu ^uj^)t~ ixt^l2\ 



dt 



\j (l^ + (c^ - c^o + Xtf) 
exp [i {uj - cjo) t' + »X*'^/2] 
y (7' + (c^ - c^o + Xt'f) 



Ca{t') 



dt'K{t,t')Ca{t'), 



(18) 



p{uj)g{uj,t)g{uj,t') = 



I" V(7^ + (^fl(^7 1) - ^0)^) 
1 

(16) 

To proceed with the analysis, a specific form for the time 
dependence of the modes must also be chosen. In this 
paper we choose a simple case where all the bath-mode 
frequencies increase linearly at the same rate, i.e. a linear 
chirp, for which 



(17) 



Some motivation for this will be given in section |Vl If 



i.e. an integro-differential equation with the kernel given 
by 



K{t,t') =^ 



exp [—1 {lu — cuq) t — 



°° y(72 + (c.-c.o+xt)') 
exp [i (w - Wo) t' + ixi'V2] 
y (7' + (c^ - C.0 + Xt'f) 



(19) 



We also note that if we change variables to t — t — t' , and 
A = cj — cjo-|-x(t + i')/2, equation (fT5|) can be expressed 
as 



9ca(t) 
dt 



dr 



dA- 



-iAt 



V[7^ + (A + xt/2)2] [72 + (A - xr/2)2] 
I 



(20) 



III. EVOLUTION OF A STATIC RESERVOIR 



a pseudomode amplitude B(t) is defined as 



We are investigating the effects of manipulating a 
reservoir of field modes which will interact with an atom. 
That is, the reservoir changes in time. However, it will be 
helpful to discuss first the behaviour of the static prob- 
lem. Here, none of the reservoir mode frequencies vary 
with time (x = 0), so equation (flS)) simphfies to 



dcgjt) 
dt 



dt' 



duj- 



-00 + {UJ — LOo)' 

M*-'')ca{t')dt' 



rCait') 



= - / Ks(t,t')Ca{t')dt' . 



(21) 



Equation (I2ip can, for example, be solved for the dynam- 
ics of the atom by means of Laplace transforms. We may 
also use the method of pseudomodes [211] . For this system 



We-^' / e^*' Ca{t')dt' , 



(22) 



and will represent the dynamics of the field. The dynam- 
ics of Ca{t) are then contained in the following two equa- 
tions which can be exactly solved and represent equation 
(HJ) [2i|: 



dCa{t) 

dt 
dB{t) 
dt 



iDB{t) 

iDca{t)-jB{t). 



(23) 
(24) 



If the frequency scale associated with the reversible dy- 
namics of the atom-field coupling D exceeds that for irre- 
versible decay 7, then there will be a resonant exchange 
of energy between atom and cavity (Rabi oscillations) as 
a photon emitted by the atom is repeatedly emitted and 
re-absorbed. This is the regime of strong coupling. We 
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note that in the weak couphng hmit the decay of the cav- 
ity field is seen to be given by 7. To improve the clarity of 
later discussions, we will now analyse the dynamics sepa- 
rately in the strong-coupling and weak-coupling regimes. 



A. Strong Coupling 

In the regime of strong coupling to the static reservoir 
(L* 7), the exact solution to equations (|23l24p can be 
written as 



Cait) 



where [22 



cos 



■r 



(25) 



(26) 



The population of the excited atomic state then under- 
goes Rabi oscillations at angular frequency (« 2D), 
under an envelope decaying at rate 7. The atomic decay 
at a rate 7 is a feature of strong coupling to the cav- 
ity, D > 7/2 in equation (j26p . and the chosen reservoir 
structure p4)) . 

For later comparison we will determine the bath spec- 
trum S{u!,t) defined in equation (|13p . In principle this 



J 



can be found by inserting the solution ([25)) into equa- 
tion ([7|) and integrating. The resulting expression for 
S{uj, t) is very complicated and here we will instead derive 
this spectrum in the strong coupling limit where D 3> 7. 
Thus we take the solution (|25p to be approximated by 



(27) 



Ca{t) « e~'^*/2cos(m/2). 



Substituting the expression into ([7|) and integrating (us- 
ing Cfc(O) = 0) gives 



Ck{t) 



gk_ / cxp [i{uik - cjQ - r^/2 + i^/2)t\ - 1 
Y \ ujk -i^Q- f2/2 + «7/2 

exp [i{u)k — Wo + ^^/2 + i-f/2)t] — 1 
Wfc - Wo + ^/2 + 27/2 



(28) 



We now substitute ([^5[) into the definition for the spec- 
trum, equation ()13p . and obtain an expression for S{uj, t) 
which simplifies in this strong coupling limit. This is 
because S{u!,t) is small away from either lo ~ loq or 
wRiwo±r2/2. As a result, if we also use the replacement 
ujk and the reservoir structure function the 

spectrum can be simplified to: 



S{co,t) 



7/2 



27r 



(uj-LUo- r!/2)2 + (7/2)' 
7/2 



27r 



(w - Wo + n/2)^ + (7/2)^ 



1 + e"''* - 2e~''*/2 cos ((w - loq - n/2) t) 
1 + e"''* - 2e"^*/^ cos ((w - wq -f VL/2.) t) 



1 



e-^* sin^ {nt/2.). 



TT (w - Wo) + 7^ 

In the long time limit the atom has completely decayed and the initial energy is distributed in the bath as 



(29) 



5(w,t ^ 00) = -!- 



7/2 



7/2 



(w - Wo - VL/2Y + (7/2)2 (cc; - wo + n/2)2 + (7/2)' 



(30) 



(See also Ref. ^23,].) 



It is seen from equation ([^5[l that the atom interacts with 
modes in the approximate range wq — 51/2 to wq -|- Sl/2 
(in the sense that only modes in this range become sig- 
nificantly changed via interaction with the atom). As 
the interaction progresses the reservoir mode occupation 
probabilities develop away from the atomic freq uency 
Wo at Wo ± ri/2 (the vacuum Rabi-splitting, [24 [25j). 
For the intermediate coupling case, where the tails of 
the Lorentzians in equation ()29p overlap, the population 
appears to move outwards (figure [2]). However, equa- 



tion ()29p is still approximately valid. For this intermedi- 
ate coupling the system does undergo Rabi oscillations, 
[D > 7/2), without the strong coupling condition £) ^ 7 
being satisfied. In this same regime we find that chirping 
the reservoir mode frequencies has interesting effects as 
we will see in the next section. This is because the occu- 
pied modes at frequencies less than wo may be brought 
back into resonance with the atom through the reservoir 
frequency manipulation and thus affect |ca(/:)l^ via equa- 
tion ©. 
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RESULTS FOR A TIME-DEPENDENT 
RESERVOIR 



[uj - wo)/7 

FIG. 2: Occupation of reservoir field-modes S{lu, t) at four 
different times. The times are chosen so that the atomic ex- 
cited state probability, |ca(t)p is either zero {t — tv/Q,, Stt/JI), 
or at a local maximum {t = An /Q,,9,n /Q). The reservoir is 
static (x = 0) and modes in the range cjq — to (^o + ^^/2 
are occupied over the four Rabi cycles shown here. The pa- 
rameter D = 87 (strong coupling). The calculation is based 
on the substitution of equation (I25p into ((T]). 



B. Weak Coupling 

For reference later we note below the behaviour of the 
atomic dynamics in the weak coupling regime {D ^7). 
The exact solution to equations (|23l24p shows damping 
at two different rates given approximately as D"^ and 
7. The damping at rate 7 only has a small effect at short 
times. An alternative approach notes that the kernel 
Ks{t,t') in equation (pij) is negligible away from t' « t. 
Then the coefficient Ca{t') can only contribute to the inte- 
gral at t' « t, so we may make the Markov approximation, 
and equation pT|) simplifies to 



dt 



Cait), 



where, 



m =^ (I - e-^o 

7 



(31) 



(32) 



This time-dependent decay rate gives atomic dynamics 
which are correct to second order in (D/7)^ and approx- 
imates the short-time dynamics as well as the long time 
behaviour. For the latter, we see that the term e~'^* in 
r(t) decreases rapidly over the decay-time of the atom, 
so we may then let 



2D^ 



(33) 



resulting in exponential time-dependence for the atom: 
Ca{t) w exp(— D^t/7). The atom and reservoir become 
entangled and the amplitude Ca{t) decays at a rate which 
is always less than the strong coupling case, equation 

(123. 



IV. 



Exact, analytic solution of equation ([T8|) has not 
proved possible. Square roots in the denominator in- 
troduce branch cuts in the complex angular frequency 
plane, and this creates difficulties in more than one ap- 
proach. Therefore the solution to this problem must be 
numerical. However, in certain limits, an approximate 
analytical solution is possible, and the results are pre- 
sented in sections HVXnilV A 31 and UVBl 

Where a numerical solution for the amplitudes Ca{t) 
and Cfc (t) has been required, it was obtained by integrat- 
ing ^ and ([7]) for a discrete micro-bath of field modes, 
with the form (|14p for the coupling constants, and linear 
chirp (|17p for the time-dependence of the reservoir mode 
frequencies. The numerical method used was a fourth- 
order Runge-Kutta stepwise integration, with adaptive 
step size. A suitable density of states was chosen and 
the results were checked to make sure the dynamics of 
Ca(i)P and |cfc(i)p were insensitive to this choice. (Typ- 
ically about ten modes within one half-width 7 of the 
reservoir structure function ([T3]) was sufficient.) Results 
of this simulation were tested in the static (x = 0) case 
against the exact solution afforded by the method of 
pseudomodes (equation pS)) and reference [2lj), and for 
very high chirp rates against analytical predictions as 
discussed later in this paper (equation (001)). 



A. Strong Coupling 

The effects of reservoir chirp are examined in three 
regimes: low, intermediate and high chirp-rates. The 
classification depends on the rate of change of the mode 
frequencies x |26[, the time-scale of the atomic dynamics, 
and also the width of the occupied mode spectrum. As 
discussed in section IIIIl even in the absence of frequency 
modulation, the time-scale of the atomic dynamics de- 
pends on ratio of the weight D to the width 7 of the cou- 
pling spectrum. For strongly coupled systems {D ^7), 
the relevant time-scale is the Rabi period 27r/J7, and since 
modes in the frequency range ujq — D,/2 to loq + in- 
teract significantly with the atom, a relevant parameter 
is the ratio of the frequency increase that each mode ex- 
periences during one Rabi period to the frequency range 
half-splitting ^l/2. Thus we introduce 



(34) 



as the parameter we will use to characterise the differ- 
ent types of behaviour. Then the specification of only 
two parameters in the model completely determines the 
dynamics: i.e. D/"f (the strength of the atom-cavity cou- 
pling) and the scaled chirp rate ^. Within this plane in 
parameter space, three categories of solution are found, 
as shown in figure [31 Here 'low' chirp corresponds to 
^ ^ 1, 'high' chirp to ^ ^ 1, and all three categories will 
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when t ^ i.e. the Lorentzians under the square-root 
separate when x(^~^0 7- Making the Markov approx- 
imation, the integro-differential equation (fT5)l becomes 



dt 



Ca{t) 



exp 



\j (7^ + {uj -ujo + xty 



X / dt' 



* ,exp[iiu;-LUo)t' +'^t'^] 



72 + (cj - + xt'Y 



Ca{t) j dt'K{t,t') 



(35) 



FIG. 3: Parameter space portrait highlighting the different 
character of solutions to equation (|18|) in different regions. 
The specific cases in figures I4I9I are indicated. The ana- 
lytic result of equation (|39|l holds in the weak-coupling regime 
{D < 7/2), and also in the 'high-chirp' regime S> 1), irre- 
spective of D/'y. 



be discussed separately below. Note that the definition 
of ^ only makes sense for strong couphng, D > 7/2, as 
this is the region in which Q is real and non-zero (see 
equation ([26]) ). However, in figure[3l to include both the 
strong and weak coupling regimes, we use irx/ D"^ for the 
scaled chirp-rate, instead of ^, above. For most of the 
plane shown, il « 2D, and so this difference is minor. 



1. High Chirp Rates 

In the limit <^ ^ 1, each mode of the reservoir is ef- 
fectively only coupled to the atom for a small fraction 
of one Rabi period. During this time, the occupation of 
a mode may increase slightly through equation ([S]), but 
does not grow large enough to re-populate the atomic 
excited state via equation ([7]). Therefore, although the 
atom-reservoir coupling is strong, Rabi oscillations are 
not observed. Manipulating the reservoir frequencies can 
change the type of atomic state behaviour from non- 
Markovian to Markovian. A specific example is shown 
in figure [H which lies in the region of high chirp rate 
indicated on figure [31 In this case there is a dramatic 
difference in the evolution of the atomic population be- 
tween the static reservoir (dashed oscillating curve) and 
the chirped reservoir (exponentially decaying curve). 

To investigate this we note that when we chirp the 
mode frequencies at a high rate, the two-time product 
p{u)g{u!,t)g{(jj,t') of equation (|16p is only non-negligible 



where we have defined a time-dependent decay rate as 



r{t) = 2 / K{t,t')dt'. 



(36) 



Time-dependent decay rates similar to r(i) also arise 
in the time-convolutionless projection operator technique 
[27t in the second order approximation. We already en- 
countered a similar time-dependent decay rate in the 
weakly coupled static reservoir in equation p2p . The 
decay rate T(t) tends to a constant value as t — > 00, 
which corresponds to a Markovian limit. This limit is 
given by 




FIG. 4: Conversion from non-Markovian to Markovian dy- 
namics by the introduction of reservoir chirp. In this exam- 
ple the dynamics were non-Markovian (dashed line) when the 
atom interacted with a static reservoir, as manifested by the 
Rabi oscillations. Chirp in the reservoir frequencies resulted 
in Markovian dynamics for the atom, and an exponential de- 
cay of |ca(t)|'^ (solid line). The parameters used were D — 



and for the chirped case, x = 4OO7 , which gives ^ 



20. 



Overlaid on the solid line (and nearly indistinguishable from 
it) , is the analytic result of equations (|39l40p (shown dotted) 
where Too ~ 37. 
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Too = lim r{t) 

t—>OG 



= lim 



2D^-f /■* p exp {i {UJ-CUO + xt' /2) t') 
exp (-i (w - Wo + Xi/2)) 



dwdt' ^ . (37) 



To evaluate this we may change to variables = x{t — t')/2 and A = u; — ujo + x{t + 1')/2, so that 
- ^dA^TdA- exp(-2.AA,/,) 



Jo J~oo y/i^^ + {A + A^r) {j^ + {A - Ad^) 
'^'^ r r _ oxp(-2.AA./x) ^33^ 



J-oo J-oo v/(7' + (A + ^<i)2) (72 + (A ~ Ad)2) 



Then changing variables again, to C/ 
y = ( A — Ad) /7, we finally obtain 



(A + Ad)/7 and 



exp 



Kq 



nx 



4X 



2 



dU 



exp 



(^) 



-17 



(39) 



Here iiTg is a modified Bessel function of zero order and 
Too is necessarily real, since Kq{z*) = Kq{z)* . Further, 
we note that Too is a monotonically decreasing function 
of X- We will later argue that this expression for Too 
is valid for ^ » 1 in the strong coupling case since the 
separation of the Lorentzians in equation (|35p . for which 
xit~t') ~ 7j must take place on a time-scale shorter than 
the time-scale for the atomic dynamics. (This would be 
1/D in the case of strong-coupling and no chirp.) With 
such high chirp, resonant modes which start to become 
occupied are moved far from resonance before they can 
act back on the atom, and are replaced by empty ones. In 
this way, the atom effectively sees an empty reservoir and 
obeys Markovian dynamics. A solution for the behaviour 
of Ca{t) follows trivially: 



c,(t) «exp(-roo t/2). 



(40) 



Figure [H] shows the agreement between this approximate 
analytic result and a numerical simulation which used a 
bath of discrete field modes and did not make the Markov 
approximation. In fact, the result (|40| also holds in the 
weak-coupling regime, independent of the rate x, as dis- 
cussed in sections flV Bl and IVl For this reason as x — ^ 
the static reservoir result Foo = 21)2/7, equation ([55]) is 
recovered. However, in the strong coupling case, this 
limit cannot be attained with equation (I34p . since, even- 
tually, the constraint ^ ^ 1 is violated. In addition, 
we know that for x = we should find the static decay 



result 7 implied by equation pSj) . Nevertheless, since 
2D'^ /j > 7 in the strong coupling regime, it is straight- 
forward to have an enhanced decay rate (i.e. Foo > 7) 
if ^ is greater than unity, but not too high. In figure [S] 
we can find this for the low chirp part of the £> = 27 
curve (12 < S 38), since only then is F00/7 > 1 

(with ^ ^ 1); for the rest of the curve the decay is inhib- 
ited relative to x = 0. However, if D/j is even modestly 
increased the region of enhanced decay is increased sub- 
stantially. It is interesting to note the equation ((39| can 
be re-written as 



F00/7 = 2iD/^f 



2\Ko{i/x)\ 



(41) 



where x — 4x/7^. In this expression the inequality 
2\Ko{i/x)\'^ /ttx < 1 holds, with equality as x (i.e. x) 
tends to zero. Since the factor D/j is always greater 
than unity for strong coupling, we may have either en- 
hanced or inhibited Markovian decay, depending on the 
value of X- Equation (pijl can also be used to motivate 
the constraint ^ » 1 introduced above. Since the time- 
scale tat for atomic dynamics is now I/Fqo, and the factor 
2\Ko{i/x)\'^ /ttx is always less than (or equal to) unity, it 
follows that tat > 1/20^. Since the time- width of the 
kernel is approximately 7/Xi the system will be Marko- 
vian if 7/x "C tat- Since ^ ^ ttx/D"^ for D > 7, this will 
always be true if ^ ^ 2it. 

For very large chirp, compared to 7^, we can expand 
the Bessel functions in equation to obtain the ap- 
proximate result 



/7->(^)^^[ln(4x/7^)]^ (42) 



which slowly decreases as x/7^ is increased. 



15 30 45 60 

Chrip rate, xjl^ 

FIG. 5: Dependence of decay rate Too on the chirp-rate 
X- The dashed lines show the prediction of equation (I39|l . 
whereas the crosses show the decay rates extracted from nu- 
merical solution (without making the Markov approximation) 
by fitting an exponential decay to the data. The parameters 
used for the three curves were D = 7/2,7, and 27. Note 
that the fit is only expected to be good for ^ 3> 1, since for 
smaller chirp-rates than this, the behaviour of \ca{t)\^ is not 
necessarily an exponential decay and may even be oscillatory. 
Thus for the highest two values of D the position of ^ = 10 is 
marked. The chirp x enhances the decay when roo/7 > 1. 

2. Intermediate Chirp Rates 

As already discussed in section UlIAl in the static case 
the excitation of the reservoir modes moves outwards 
from ojo over several Rabi cycles towards final peaks at 
tJo i Now in the intermediate-chirp case, the fre- 

quency of the modes is increased by approximately the 
same amount over one Rabi-cycle. That is, the occupied 
modes are brought back into resonance with the atom 
on every Rabi-cycle. This can have a profound effect 
on the dynamics of the atomic excited state probability, 
Ca(t)P, since occupied modes with angular frequency in 
the vicinity of ujq have the strongest interaction with the 
atom, as can be seen from equation 

In this regime of intermediate chirp rate 1) and 
strong coupHng (£) ^ 7), we must rely on numerical re- 
sults. As illustrated in figure [31 two free parameters, i.e. 
D/7 and the rate of redistribution of reservoir mode fre- 
quencies, label the dynamics of the system. A numerical 
search was conducted over this plane, and it was found 
that Rabi oscillations were inhibited over most of this re- 
gion. The explanation for this is that occupied modes of 
the reservoir are moved from resonance before they have 
a chance to re-populate the atom. A typical example is 
shown in figure [S] Additionally, for -D/7 < 15, there is 
a particular chirp-rate for which a 'recycling condition' 
is fulfilled and the vacuum Rabi-splitting in the reser- 
voir doesn't occur. A significant population may become 
trapped in a cycle of resonant exchange between atom 
and reservoir, and Rabi oscillations may persist for much 
longer than in the static case. The atomic behaviour 
is shown in figure [T] We also note that the oscillation 
frequency increases with increasing chirp-rate as is seen 
empirically from simulations. In figure [5] we see the cor- 
responding spectrum of excitation in the reservoir for two 



1 2 3 4 5 

It 

FIG. 6: Figure showing enhanced decay of the envelope of 
|ca(t)P caused by manipulation of the reservoir modes to 
which the atom is coupled (by a linear chirp). We note that 
the Rabi period is modified by manipulation of the reservoir 
mode frequencies. Parameters are D = 87 and for the solid 
line, X = 207^. Thus ^ ~ 1.0. For comparison, the dashed 
line shows |ca(t)P for the case with no chirp {x~Q)- 
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It 



FIG. 7: Inhibition of decay and trapping of population with 
a chirped reservoir. The parameters are D = 87 and for the 
solid line, x ~ 8.47^ (which implies ^ ~ 0.42). For compari- 
son, the dashed line shows |ca(f)P for the case with no chirp 
(X = 0). 



pairs of times. In each case the structured peak at higher 
frequencies is sufficiently detuned for this reservoir pop- 
ulation to be effectively decoupled from the atomic dy- 
namics. Therefore, excitation simply moves to higher fre- 
quencies at later times, as may be seen by comparing the 
spectra after 10 and 20 Rabi cycles. This non-interacting 
peak has an area of approximately a half, which rep- 
resents the population that was 'pushed' to frequencies 
above ujq during the Rabi splitting at early times. We see 
that the left hand peak, which is centred at the atomic 
transition, remains on resonance as individual mode fre- 
quencies pass through this region. It does invert peri- 
odically during the time evolution, but this population 
exchange, or recycling, between the atom and locally res- 
onant reservoir modes is very stable. For example, we 
can see that at the same part of the cycle the shape of 
this peak is relatively unchanged between (a) and (c) [and 
between (b) and (d)] in figure [5] This behaviour is quite 
different to that seen for the static reservoir in figure [21 
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3. Low Chirp Rates 



FIG. 8: The reservoir excitation spectrum during the sta- 
ble oscillations of figure [7] The horizontal axis shows the 
scaled time-dependent detuning from the atomic transition 
frequency (vertical dashed line), where u>r = uj + xt- For 
(a-d) "ft ~ 3.81, 4.00, 7.60, 7.79, and other parameters are 
as for figure [T] i.e. D — 87, x ~ 8.47'^. The times selected 
correspond to approximately: (a) 10; (b) 10.5; (c) 20; and (d) 
20.5 Rabi oscillations. 



In the regime of strong coupling and very low chirp 
^ 1) the basic behaviour of the system will be to 
perform Rabi oscillations of the type found in equation 
([25|) . The weak chirp will be expected to perturb the 
dynamics only slightly. In the example shown in figure 
[9] the difference between the chirped and unchirped dy- 
namics is very small, although it may be noticed the the 
chirped bath results in a slightly raised Rabi frequency. 
This frequency shift might be observable if the decay rate 
is low enough that a large phase shift accumulates in the 
Rabi oscillations. To evaluate the size of this effect we 
examine a perturbative solution found by expanding the 
denominator of equation (|20p in powers of x- Thus, for 
the kernel of equation ((20|) 



Kir) 



TT 
12, 



dA 



dA 



exp(— lAr) 
Vb^ + (A + Xt/2)2] [72 + (A - 



72 + A2 



4 (72-)-A2)2 



+ 



(43) 



to order x^- The method of pseudomodes (see section ITTTl 
and Ref. ,21]) could be used to analyse this kind of kernel 
structure, but here we will simply examine the Laplace 
transform of equation (I^O)) to find the perturbation of the 
Rabi frequency. To do this we first perform the integral 
over the detuning A in equation so that with the 
approximate kernel (and given r > 0) 



If we momentarily set x to zero we see that, if we 
substitute equation into equation ([l5|) . the poles 
in the right hand side of equation (|45l) are located at 
s = —7/2 ± iil/2, where ft is given in equation (|26p . 
This represents the damped oscillations of equation (|25p . 
If we now admit a finite value of x, we seek the perturbed 
locations of the poles in the approximate form 



K{t) = D'^e-^'^ 



2 2 



(44) 



If we Laplace transform equation ([20]) and solve for Ca{s), 
the transform of Ca{t), we find 



Ca(s) = 



Kis)' 



(45) 



where K{s) is the transform of equation ((43l) and we let 
Ca(0) = 1. Using the approximate kernel (j44|) wc will find 
that 



Kis) 



s + 7 



1 - 



X 



872(5 + 7)2 
3x' 3x2 



87(5 + 7)3 2(s + 7)i 



(46) 



7 , .^^ 
-2^*2 



(47) 



The factor S is found by substitution into the equation 
s + K{s) — 0. Since the modification of the decay rate 
is small, and not of interest, we examine the imaginary 
part of S to find that Im((5) ~ ±l/{8il'j'^) and then the 
new Rabi frequency fi' is approximately 



r?' - 1 



X" 



4f^2-^2 



(48) 



In the above we also used the strong coupling approxima- 
tion (£',n) ^ 7. Even if the correction to the Rabi fre- 
quency is small, the effect may be noticeable since there 
will be an accumulated phase shift over time. Since the 
largest usable time in this limit is of order I/7, it follows 
that the largest phase difference is approximately given 
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by 



(49) 



Thus even if the first term in brackets is small, since it 
is the correction to the Rabi frequency in equation (^5]) . 
the strong coupling regime ensures the second term is 
large and so the resulting phase difference may be of a 
reasonable size. Note that the effect will be made more 
visible, not only by increasing the chirp but also by 
reducing the strong atom-cavity coupling. 




FIG. 9: For low chirp-rates, the atom experiences coupling 
to the same set of modes for several Rabi cycles and, conse- 
quently, the effects of the chirp are small. Parameters chosen 
were D = 87, and for the solid line, x = 27^, which corre- 
sponds to ^ ~ 0.1. For comparison, the dashed line shows 
|ca(t)P for the case with no chirp (x = 0). 



B. Weak Coupling 

In the limit D ^ 7, the behaviour of |ca(i)P is Marko- 
vian. This remains true, irrespective of how the reser- 
voir mode frequencies are manipulated. Therefore, for 
weakly-coupled reservoirs, the analysis of section HV A 11 
culminating in equation (|40l) . is valid for all chirp-rates. 
The region of weak coupling corresponds to the horizon- 
tal strip at the bottom of figure [31 Note that the expan- 
sion [281 



1 



hm {Ko{z)}^J-( 



1 

8z 



128z2 

(50) 

implies that the decay constant Too tends to the result of 
the static reservoir, equation (j33p . in the low-chirp limit, 
as it should. We can also find the leading correction, so 
that as X — * 



2D^ 



1 



2x' 



(51) 



V. CONCLUSION 

We have investigated the effects on the dynamics of 
an open quantum system of a particular manipulation of 



the reservoir mode frequencies. While other authors have 
discussed reservoirs with time-dependent properties, we 
have here considered the subtle case where all macro- 
scopic properties of the reservoir remain static, while 
the microscopic structure of the reservoir changes. This 
problem is intriguing, since the effects are non-trivial, 
and cannot easily be guessed in advance. It is tempting 
to think that since the only difference with the usual, 
static case is a redistribution of reservoir mode popu- 
lations, only non-Markovian systems would be affected. 
This is not the case, however, and an effect can be seen 
for both Markovian and non-Markovian systems. 

We have characterised the nature of the dynamics in 
three different regions of parameter space (see figure [Hj) 
and found conditions for each type of behaviour in section 
IIVI It is shown that manipulation of the reservoir mode 
frequencies may alter the nature of the dynamics (from 
non-Markovian to Markovian behaviour) when the chirp- 
rate is high (^ ^ 1). Analytical results are presented for 
the decay in this case (equations and (dOj))- The 
resulting decay may be enhanced or inhibited, compared 
to the un-chirped case. In the case of strong coupling 
and low chirp, the chirp manifests itself as a frequency 
shift in the Rabi oscillations (equation ([48|)). 

In order to realise the effects of section IIVI in the con- 
text of optical cavity QED it would be necessary to con- 
struct a high-Q cavity and a reservoir whose mode fre- 
quencies could be manipulated as described above, i.e. 
without altering the envelope of the coupling constants 
near the atomic resonance. Ideally, if the resonant fre- 
quency of the cavity is to remain unchanged (at a;o/27r), 
its cavity mirrors must be fixed, too. Since there will be 
losses through the cavity mirrors, i.e. a coupling to the 
external environment, the bath mode frequencies may be 
altered by changing the properties of that external en- 
vironment. This would be sufficient to isolate the basic 
effect from, for example, changing the position of the cav- 
ity resonance (as in Ref. fl3]). The most simple way of 
doing this would be to have an inner and outer cavity, 
with lengths i and L. We can suppose that the outer mir- 
ror can be moved by a piezoelectric actuator, or similar 
device, which will change the length L. Then for small 
fractional changes in the outer cavity length we can write 



'dt 



Lo dL 
'L~dt 



luq dL 

'T~dt' 



so that 



X ~ — r 



(jjQ dL 
T~dt' 



(52) 



(53) 



To illustrate the potential size of any observable ef- 
fects, we consider some specific cases. First we take 
an optical atomic transition with angular frequency 
loq/2tt ^ 3.5 X 10^^ Hz and a 40 /im inner cavity for 
which 7/27r ~ 4.1 MHz and i:»/27r - 34 MHz [H HI]. 
This is in the strong coupling regime. For the outer cavity 
we let L = 1cm, and we assume here that the inner cav- 
ity is diffractively coupled to the outer cavity to enhance 
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the mode density, whilst maintaining a strong directional 
property. If we can move the out-most mirror at a speed 
of 0.1 ms^^ for the short time necessary, we find that 
X ^ 2.2 X 10^^ s"^. Then ^ « 1.52, which is enough to 
enter the intermediate chirp region indicated in figure [3] 
and discussed in section llV A2I If we would increase the 
mirror speed to 0.65 ms^^, ^ increases to approximately 
ten, which places us in the high chirp, strong coupling 
regime of figure [3l In this case equation ([39)l applies and 
we obtain an enhancement of decay: Too/^ = 5.05. This 
enhanced decay can be increased even further if, for ex- 
ample, the cavity decay rate is reduced by a factor of 
ten: in that case Too/"/ = 13.6. As an example of in- 
hibited decay we can simply reduce D so that D ^ ■y. 
In this case we are in the weak coupling regime and the 
inhibition of decay relative to the static reservoir case is 
independent of D and given by a factor of approximately 



(2i?2/7) 



9.92 X 10"^ 



(54) 



when 7/27r ^ 0.41 MHz and the speed is 0.65 ms ^, as 



above. 

We have seen that many of the effects described here 
could be tested with currently available experimental pa- 
rameters. For an atom weakly coupled to a high-Q cavity, 
the inhibition of decay and decoherence is considerable, 
with the possibility of suppressing the decay rate by a 
factor of 10"^. In the case of strong coupling, atomic de- 
cay can be inhibited (for high chirp x) or enhanced. The 
examples given show an enhancement of decay by a factor 
of 5 or 14. This factor can be increased if the ratio D/y is 
increased along with the smallest cavity chirp consistent 
with ^ 1. An attractive feature of this model is that 
the suppression, or enhancement, could be turned on and 
off at will by moving or not moving the mirror of an outer 
cavity. In the case of enhanced decay, this means that the 
controlled motion of an outer cavity mirror can extract 
energy from the inner cavity-atom system. 

This work was supported in part by the UK Engineer- 
ing and Physical Sciences Research Council. 
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